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^ ' In these lectures, a simple introduction to the phenomenon of diffraction 

, in deep inelastic scattering and its theoretical description is given. While 

' the main focus is on the diffractive structure function , some issues in 

' diffractive vector meson production are also discussed. 

a^ ■ 
o 

^ ' 1. Introduction 

o> ■ 

r~| ' In the following, the term 'diffractive electroproduction' or 'diffractive 

Qh! deep inelastic scattering (DIS)' is used to characterize those processes in 

Q^' small-x DIS where, in spite of the high virtuality of the exchanged photon, 

O . the proton target remains intact or almost intact. Thus, the term is used 

synonymously with the expression 'rapidity gap events in DIS'. 

Historically, the term diffraction is derived from optics, where it de- 
^ I scribes the deflection of a beam of light and its decomposition into com- 

^ • ponents with different frequencies. In high energy physics, it was originally 

used for small-angle elastic scattering of hadrons (Fig. la). If one of the 
hadrons, say the projectile, is transformed into a set of two or more fi- 
nal state particles, the process is called diffractive dissociation or inelastic 
diffraction. Good and Walker have pointed out that a particularly intuitive 
physical picture of such processes emerges if the projectile is described as 
a superposition A + B oi different components which scatter elastically off 
the target [1] (Fig. lb). Since the corresponding elastic amplitude is differ- 
ent for each component, the outgoing beam will contain a new superposition 
aA + PB of these components and therefore, in general, new physical states. 
These are the dissociation products of the projectile. 

Even at very high energy, the above processes are generically soft, i.e., 
the momentum transfer is small and the dissociation products have small 
p±. Therefore, no immediate relation to perturbative QCD is apparent. 
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Figure 1: Elastic scattering (a) and difFractive dissociation (b) in the phys- 
ical picture of Good and Walker [1]. 

In contrast to these soft processes, diffractive DIS is an example of hard 
diffraction, the hard scale being the virtuality of the exchanged photon. 
However, this is not the first hard diffractive process that was observed. 
Hard diffraction was first investigated in diffractive jet production at the 
CERN SppS collider in proton-antiproton collisions [2] . Although one of the 
hadrons escapes essentially unscathed, a high-pj^ jet pair, which is necessar- 
ily associated with a high virtuality in the intermediate states, is produced 
in the central rapidity range (Fig. 2). The cross section of the process is 
parametrically unsuppressed relative to non-diffractive jet production. This 
seems to contradict a naive partonic picture since the colour neutrality of 
the projectile is destroyed if one parton is removed to participate in the 
hard scattering. The interplay of soft and hard physics necessary to explain 
the effect provides one of the main motivations for the study of these 'hard 
diffractive' processes. 




X 



± 



Figure 2: Diffractive dissociation in hadron-hadron collisions. If the diffrac- 
tive final state X contains two high-px jets, the process is called diffractive 
jet production. 

Here, the focus is on diffractive electroproduction, which is another ex- 
ample of a hard difFractive process. This process became experimentally 
viable with the advent of the electron-proton collider HERA (cf. [3]), where 
DIS at very small values of the Bjorken variable x can be studied. In the 
small-x or high-energy region, a significant fraction of the observed DIS 
events have a large rapidity gap between the photon and the proton frag- 
mentation region [4,5]. In contrast to the standard DIS process X 
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(Fig. 3a), the relevant reaction reads 7*p XY (Fig. 3b), where X is a 
high-mass hadronic state and Y is the elastically scattered proton or a low- 
mass excitation of it. Again, these events are incompatible with the naive 
picture of a partonic target and corresponding simple ideas about the colour 
flow. Naively, the parton struck by the virtual photon destroys the colour 
neutrality of the proton, a colour string forms between struck quark and 
proton remnant, and hadronic activity is expected throughout the detector. 
Nevertheless, the observed diffr active cross section is not power suppressed 
at high virtualities with respect to standard DIS. 
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Figure 3: Inclusive (a) and diffractive (b) electroproduction. 

These notes are organized as follows. In Sect. 2, the basic experimental 
observations and the diffractive structure function are discussed. Sec- 
tion 3 describes the theory of ^2^, emphasizing the semiclassical approach 
and its relation to the concept of diffractive parton distributions. In Sect. 4, 
some topics in diffractive meson production, in particular the factorization 
of the hard amplitude, arc described. A brief summary is given in Sect. 5. 

The reader may consult [6] for a more detailed review. 

2. Basic observations and the diffractive structure function 

In this section, the kinematics of diffractive electroproduction at small 
X is explained in some detail, and the notation conventionally used for the 
description of this phenomenon is introduced. The main experimental obser- 
vations are discussed, and the concept of the diffractive structure function, 
which is widely used in analyses of inclusive diffraction, is explained. 

2.1. Kinematics 

To begin, recall the conventional variables for the description of DIS. 
An electron with momentum k collides with a proton with momentum P. 
In neutral current processes, a photon with momentum q and virtuality 
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Figure 4: Diffractive electroproduction. The full hadronic final state with 
invariant mass W contains the elastically scattered proton and the diffrac- 
tive state with invariant mass M. 



q'^ = —Q^ is exchanged, and the outgoing electron has momentum k' = k — q. 
In inclusive DIS, no questions are asked about the hadronic final state Xw, 
which is only known to have an invariant mass square W'^ = {P + q)^. 
The Bjorken variable x = Q'^ /{Q'^ + W'^) characterizes, in the nai'vc parton 
model, the momentum fraction of the incoming proton carried by the quark 
that is struck by the virtual photon. If a; <C 1, which is the relevant region in 
the present context, Q is much smaller than the photon energy in the target 
rest frame. In this sense, the photon is almost real even though ^ 
(where A is some soft hadronic scale) . It is then convenient to think in terms 
of a high-energy j*p collision with centre-of-mass energy W. 

Loosely speaking, diffraction is the subset of DIS characterized by a 
quasi-elastic interaction between virtual photon and proton. A particularly 
simple definition of diffraction is obtained by demanding that, in the 7*p 
collision, the proton is scattered elastically. Thus, in diffractive events, the 
final state contains the scattered proton with momentum P' and a diffractive 
hadronic state Xm with mass M (see Fig. 4). Since diffractive events form 
a subset of DIS events, the total invariant mass of the outgoing proton and 
the diffractive state Xm is given by the standard DIS variable W. 

The following parallel description of inclusive and diffractive DIS sug- 
gests itself. 

In the former, virtual photon and proton collide to form a hadronic state 
Xw with mass W. The process can be characterized by the virtuality 
and the scaling variable x = Q"^ /{Q^ + W^), the momentum fraction of the 
struck quark in the naive parton model. 

In the latter, a colour neutral cluster of partons is stripped off the proton. 
The virtual photon forms, together with this cluster, a hadronic state Xm 
with mass M. The process can be characterized by Q^, as above, and by 
a new scaling variable P = Q^/iQ"^ + M^), the momentum fraction of this 
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cluster carried by the struck quark. 

Since diffraction is a subprocess of inclusive DIS, the struck quark from 
the colour neutral cluster also carries a fraction x of the proton momentum. 
Therefore, the ratio = x/(3 characterizes the momentum fraction that the 
proton loses to the colour neutral exchange typical of an elastic reaction. 
This exchanged colour neutral cluster loses a momentum fraction (3 to the 
struck quark that absorbs the virtual photon. As expected, the product 
X = is the fraction of the original proton's momentum carried by this 
struck quark. Since the name pomeron is frequently applied to whichever 
exchange with vacuum qiiantum numbers dominates the high-energy limit, 
many authors use the notation xjp = ^, thus implying that the proton loses 
a momentum fraction ^ to the exchanged pomeron. 

Therefore, x, and f3 or, alternatively, x, and ^ are the main 
kinematic variables characterizing diffractivc DIS. A further variable, t = 
{P — P')^, is necessary if the transverse momenta of the outgoing proton 
and the state Xm relative to the j*P axis are measured. Since the proton 
is a soft hadronic state, the value of \t\ is small in most events. The small 
momentum transferred by the proton also implies that M <C W. 

To see this in more detail, introduce light-cone co-ordinates q± = qo± qs 
and q± = {qi-q2)- It is convenient to work in a frame where the transverse 
momenta of the incoming particles vanish, q± = P± = 0. Let A be the 
momentum transferred by the proton, A = P — P', and nip = P^ = P'^ the 
proton mass squared. For forward scattering, = 0, the relation 

t = = A+A_ = -fml (2.1) 

holds. Since ^ = {Q^ + M"^) / {Q^ + W^) , this means that small M implies 
small |t| and vice versa. Note, however, that the value of \t\ is larger for 
non-forward processes, where t = A_|_A_ — A^. 

So far, diffractive events have been characterized as those DIS events 
which contain an elastically scattered proton in their hadronic final state. 
An even more striking feature is the large gap of hadronic activity seen in the 
detector between the scattered proton and the diffractive state Xm- It will 
now be demonstrated that this feature, responsible for the alternative name 
'rapidity gap events', is a direct consequence of the relevant kinematics. 

Recall the definition of the rapidity y of a particle with momentum k, 

y=lln^ = lln^°±^. (2.2) 

This is a convenient quantity for the description of high-energy collisions 
along the z-axis. Massless particles moving along this axis have rapidity — oo 
or -|-oo, while all other particles are characterized by some finite intermediate 
value of y. 
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In the centre-of-mass frame of the 7*p coUision, with the z-axis pointing 
in the proton beam direction, the rapidity of the incoming proton is given 
by Up = ln(P+/mp). At small ^, the rapidity of the scattered proton is 
approximately the same. This is to be compared with the highest rapidity 
ymax of any of the particles in the diffractive state Xm ■ Since the total plus 
component of the 4-momentum of Xm is given by — x)P+, and the pion, 
with mass m,.,^, is the lightest hadron, none of the particles in Xm can have 
a rapidity above Umax = — x)P+/m7r). Thus, a rapidity gap of size 

Ay = ln(m^/(^ — x)mp) exists between the outgoing proton and the state 
Xm- For typical values of ^ 10"'^ the size of this gap can be considerable. 

Note, however, that the term 'rapidity gap events' was coined to de- 
scribe the appearance of diffractive events in the HERA frame, i.e., a frame 
defined by the electron-proton collision axis. The rapidity in this frame is, 
in general, different from the photon-proton frame rapidity discussed above. 
Nevertheless, the existence of a gap surrounding the outgoing proton in the 
7*p frame clearly implies the existence of a similar gap in the ep frame. The 
exact size of the ep-frame rapidity gap follows from the specific event kine- 
matics. The main conclusion so far is the kinematic separation of outgoing 
proton and diffractive state Xm in diffractive events with small ^. 

Without losing any of the qualitative results, the requirement of a final 
state proton P' can be replaced by the requirement of a low-mass hadronic 
state Y, well separated from the diffractive state Xm- In this case, the ar- 
gument connecting elastically scattered proton and rapidity gap has to be 
reversed: the existence of a gap between Xm and Y becomes the distinctive 
feature of diffraction and, under certain kinematic conditions, the interpre- 
tation of Y as an excitation of the incoming proton, which is now almost 
elastically scattered, follows. 

2.2. The main experimental observations 

Rapidity gaps are expected even if in all DIS events a quark is knocked 
out of the proton leaving a coloured remnant. The reason for this is the 
statistical distribution of the produced hadrons, which results in a small yet 
finite probability for final states with little activity in any specified detector 
region. However, the observations described below are clearly inconsistent 
with this explanation of rapidity gap events. 

More than 5% of DIS events at HERA were found to possess a rapidity 
gap [4,5]. The analyses are based on the pseudo-rapidity rj = — lntan(9/2), 
where 9 is the angle of an outgoing particle relative to the beam axis. Pseudo- 
rapidity and rapidity are identical for massless particles; the difference be- 
tween these two quantities is immaterial for the qualitative discussion below. 

In the ZEUS analysis, a rapidity ?7niax was defined as the maximum 
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rapidity of a calorimeter cluster in an event. A cluster was defined as an 
isolated set of adjacent cells with summed energy higher than 400 MeV. 
The measured r/max distribution is shown in Fig. 5. (Note that the smallest 
detector angle corresponds to r^max = 4.3; larger values are an artifact of 
the clustering algorithm.) 
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Figure 5: Distribution of r/max) the maximum rapidity of a calorimeter 
cluster in an event, measured at HERA (figure from [4]). 

To appreciate the striking qualitative signal of diffraction at HERA, the 
measured r/max distribution has to be compared with naive expectations 
based on a purely partonic picture of the proton. This is best done using a 
parton-model-based Monte Carlo event generator. The corresponding r/max 
distribution, which is also shown in Fig. 5, is strongly suppressed at small 
^max- This qualitative behaviour is expected since the Monte Carlo (for more 
details see [4] and refs. therein) starts from a partonic proton, calculates 
the hard process and the perturbative evolution of the QCD cascade, and 
finally models the hadronization using the Lund string model (see, e.g., [7]). 
According to the Lund model, the colour string, which connects all final state 
partons and the coloured proton remnant, breaks up via qq pair creation, 
thus producing the observed mesons. The rapidities of these particles follow 
a Poisson distribution, resulting in an exponential suppression of large gaps. 

It should be clear from the above discussion that this result is rather gen- 
eral and does not depend on the details of the Monte Carlo (note, however, 
the Monte Carlo based approach of [8], which allows for colour reconnec- 
tion). QCD radiation tends to fill the rapidity range between the initially 
struck quark and the coloured proton remnant with partons. A colour string 
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connecting these partons is formed, and it is highly unhkely that a large gap 
emerges in the final state after the break-up of this string. 

However, the data shows a very different behaviour. The expected ex- 
ponential decrease of the event number with r^max is observed only above 
^max — 1-5; below this value a large plateau is seen. Thus, the naive par- 
tonic description of DIS misses an essential qualitative feature of the data, 
namely, the existence of non-suppressed large rapidity gap events. 

To give a more specific discussion of the diffractive event sample, it is 
necessary to define which events are to be called diffractive or rapidity gap 
events. It is clear from Fig. 5 that, on a qualitative level, this can be achieved 
by an ?7max cut separating the events of the plateau. The resulting qualitative 
features, observed both by the ZEUS [4] and HI collaborations [5], are the 
following. 

There exists a large rapidity interval where the r/max distribution is flat. 
For high 7*p energies W, the ratio of diffractive events to all DIS events is 
approximately independent of W. The dependence of this ratio is also 
weak, suggesting a leading-twist contribution of diffraction to DIS. Further- 
more, the diffractive mass spectrum is consistent with a 1 /M^ distribution. 

A number of additional remarks are in order. Note first that the ob- 
servation of a flat rymax distribution and of a spectrum are interde- 
pendent as long as masses and transverse momenta of final state particles 
are much smaller than M^. To see this, observe that the plus component 
of the most forward particle momentum and the minus component of the 
most backward particle momentum are largely responsible for the total in- 
variant mass of the diffractive final state. This gives rise to the relation 
dM^/M^ = ci In ~ (i?7max; which is equivalent to the desired result. 

A significant contribution from exclusive vector meson production, e.g., 
the process j*p ^ pp, is present in the rapidity gap event sample. A more 
detailed discussion of corresponding cross sections and of relevant theoretical 
considerations is given in Sect. 4. 

2.3. Diffractive structure function 

The diffractive structure function, introduced in [9] and first measured 
by the HI collaboration [10], is a powerful concept for the analysis of data 
on diffractive DIS. 

Recall the relevant formulae for inclusive DIS. The cross section for the 
process ep eX can be calculated if the hadronic tensor, 

W^.{P,q) = < X X\MO)\P > {27r)'S\q + P-px), 

(2.3) 
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is known. Here j is the electromagnetic current, and the sum is over all 
hadronic final states X. Because of current conservation, q-W = W ■ q = 0, 
the tensor can be decomposed according to 

W^,{P,q)=(^g^,-^^ Wi{x,Q^)+(^P^ + ^q^^(^P, + ^q,^ W2{x,Q^) . 

(2.4) 

The data is conveniently analysed in terms of the two structure functions 

F2(x,Q2) = {p.q)W2{x,Q^) (2.5) 
Fl{x,Q^) = {P-q)W2{x,Q^)-2xWi{x,Q^). (2.6) 

Introducing the ratio R = Fl/{F2 — Fl), the electron-proton cross section 
can be written as 

where y = Q^/sx, and s is the electron-proton centre-of-mass energy squared. 
In the naive parton model or at leading order in in QCD, the longitudinal 
structure function Fl{x,Q'^) vanishes, and R = 0. Since R corresponds to 
the ratio of longitudinal and transverse virtual photon cross sections, ai/o-x, 
it is always positive, and the corrections associated with a non-zero R are 
small at low values of y. 

In diffraction, the two additional kinematic variables ^ and t are present. 
However, no additional independent 4-vector is introduced as long as the 
measurement is inclusive with respect to the azimuthal angle of the scat- 
tered proton. Therefore, the decomposition in Eq. (2.4) remains valid, and 
the two diffractive structure functions F^j^^^ {x, Q^,^, t) can be defined. The 
diffractive cross section reads 

dxdQ'^d^dt xQ4 y ^^2[l + i?o(4)(x,Q2,e,«)]J ' ^^'^'^'^^ 

(2.8) 

where R^ = F^ /{F2 — F['). In view of the limited precision of the data, 
the dominance of the small-y region, and the theoretical expectation of the 
smallncss of FJ^ , the corrections associated with a non-zero value of R^ are 
neglected in the following. 

A more inclusive and experimentally more easily accessible quantity can 
be defined by performing the t integration, 

F2''^''\x,Q^0 = J dtF^^^\x,Q^^,t). (2.9) 
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The main qualitative features of diffractive electroproduction, already dis- 
cussed in the previous section, become particularly apparent if the functional 
form of F^^^^ is considered (cf. the data in Figs. 13 and 14 in Sect. 3). The 

/? and dependence of F^^^^ is relatively flat. This corresponds to the 
observations discussed earlier that diffraction is a leading twist effect and 
that the mass distribution is consistent with a spectrum. The en- 

ergy dependence of diffraction is such that $,^2'^^^ (Cj P, Q"^) slowly grows as 
^ ^ 0. This growth and its possible relation to the growth of inclusive struc- 
ture functions at small x and to the growth of clastic cross sections at high 
energy is one of the most interesting aspects of diffraction (see also [11]). 

3. Theoretical approaches to 

To understand qualitatively how leading twist diffraction in DIS comes 
about, it is simplest to work in the target rest frame. Diffraction means that 
a hadronic fluctuation of the energetic virtual photon scatters off the target 
proton without exchanging colour. One can then expect the photon and the 
proton to fragment independently, leading to a rapidity gap event. 

A particularly simple physical picture is provided by the aligned jet 
model of Bjorkcn and Kogut [12]. In terms of the QCD degrees of freedom, 
a closely related formulation of diffraction was given with the two-gluon 
exchange calculations of Nikolaev and Zakharov [13]. Since the t channel 
colour singlet exchange is not hard in the bulk of the cross section underlying 
F2 , the exchange of more than two gluons is not suppressed. This can be 
systematically treated in the semiclassical approach [14,15]. 

In the Breit frame, where the target proton is fast, diffraction was his- 
torically described as DIS off a pomeron [16, 17]. A more general, QCD 
based concept treating diffraction from a Breit frame point of view are the 
diffractive parton distributions [18, 19]. As will be discussed in more detail 
below, the semiclassical approach provides a convenient framework in which 
one can intuitively understand how diffractive parton distributions emerge 
in a target rest calculation. 

S.l. Aligned jet model 

A very simple argument why, even at very high photon virtualities, 
diffractive DIS is largely a soft process was presented by Bjorken and Kogut 
in the framework of their aligned jet model [12]. 

The underlying physical picture is based on vector meson dominance 
ideas. The incoming photon fluctuates into a hadronic state with mass M, 
which then collides with the target (see Fig. 6a). The corresponding cross 
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section for transverse photon polarization is estimated by 

daT dP{M^ 



dM2 dM2 



• <t{M') , (3.1) 



where the probabihty for the photon to develop a fluctuation with mass M 
is given by 

and a"(M^) is the cross section for this fluctuation to scatter off the target. 
The above expression for dP{M'^) is most easily motivated in the framework 
of old-fashioned perturbation theory, where the energy denominator of the 
amplitude is proportional to the off-shellncss of the hadronic fluctuation, 
+ M^. If this is the only source for a dependence, the numerator 
factor is necessary to obtain a dimensionless expression. 





a) P 



Figure 6: Vector meson dominance inspired picture of inclusive (a) and 
diffractive (b) electroproduction. 

Bjorken and Kogut assume that, for large M^, the intermediate hadronic 
state typically contains two jets and that cr(M^) is suppressed for config- 
urations with high px_ (the latter effect being now known under the name 
of colour transparency). Consider hadronic fluctuations with a certain M^, 
which, in their respective rest frames, are realized by two back-to-back jets. 
Under the assumption that the probability distribution of the direction of 
the jet axis is isotropic, simple geometry implies that aligned configurations, 
defined by p]_ < (where is a soft hadronic scale) , are suppressed by 
A^/M^. If only such configurations are absorbed with a large, hadronic cross 
section, the relations (t(M^) ~ and 

dM^ ^ (M2 + Q2)2 ^•'••'^ 

follow. Thus, the above cross section can be interpreted as the total high- 
energy cross section of target proton and aligned jet fluctuation of the pho- 
ton, i.e., of two soft hadronic objects. Therefore, a similar elastic cross sec- 
tion is expected, ~ or (cf. Fig. 6b). The resulting diffractive structure 
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function reads 



^''(e,/9,Q')-^/3/e• 



(3.4) 



It is interesting that the very simple arguments outhned above cap- 
ture two important features of the HERA data: the leading-twist nature of 
diffraction and the approximate 1/^ behaviour of F2'. 



The simplest way to formulate the above intviitivc picture, where a had- 
ronic fluctuation of the photon scatters off the proton, in perturbative QCD 
is via two-gluon exchange. At leading order, the photon fluctuates in a qq 
pair and two gluons forming a colour singlet are exchanged in the t channel 
(see Fig. 7). A corresponding calculation was first performed by Nikolaev 
and Zakharov [13], where the probability for the photon to fluctuate in a qq 
pair was described by the square of the light-cone wave function of the pho- 
ton. The scattering of the qq pair off the proton was parametrized using the 
dipole cross section a{p), where p is the transverse distance between quark 
and antiquark when they hit the target. Here, the corresponding calcula- 
tions will not be described (see, however, the recent review of diffraction [20] 
emphasizing two-gluon exchange) . Note that two-gluon exchange results can 
be recovered from the semiclassical calculation described below if a Taylor 
expansion in the external colour field is performed. 



Figure 7: Diffractive electroproduction in the two-gluon exchange model. 

Gluon radiation, i.e., the fluctuation of the incoming photon in a qqg 
state was considered in [21]. 

For further interesting work related to the two-gluon exchange approxi- 
mation for see, e.g., [22-24] (cf. [25] for new developments in the BFKL 
resummation method). Note also the discussion of diffractive processes in 
the colour-dipole approach [26] and the combined analysis of diffractive and 
inclusive DIS in this framework [27]. 

The main shortcoming of the two-gluon approach is the lacking jus- 
tification of perturbation theory. As should be clear from the aligned jet 
model and as will be discussed in a more technical way below, the diffrac- 
tive kinematics is such that the t channel colour singlet exchange does not 



3.2. Two-gluon exchange 
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feel the hard scale of the initial photon. Thus, more than two gluons can be 
exchanged without suppression by powers of ag. 

Note that this is different in certain more exclusive processes, such as 
vector meson production, where it has been shown that the dynamics of the 
t channel exchange is governed by a hard, perturbative scale (cf. Sect. 4). 

3.3. Semiclassical approach 

In this approach, the interaction with the target is modelled as the 
scattering off a superposition of soft target colour fields, which, in the high- 
energy limit, can be calculated in the eikonal approximation [28]. Diffraction 
occurs if both the target and the partonic fluctuation of the photon remain 
in a colour singlet state. 

The amplitude for an energetic parton to scatter off a given colour field 
configuration is a fundamental building block in the semiclassical approach. 
The essential assumptions are the softness and localization of the colour 
field and the very large energy of the scattered parton. 

The relevant physical situation is depicted in Fig. 8, where the blob 
symbolizes the target colour field configuration. Consider first the case of a 
scalar quark that is minimally coupled to the gauge field via the Lagrangian 

C,^^l^ = {D^^Y{D^'^) with D^, = d^, + igA^. (3.5) 

In the high-energy limit, where the plus component of the quark momentum 
becomes large, the amplitude of Fig. 8 then reads 

. (3.6) 

It is normalized as is conventional for scattering processes off a fixed target. 
The expression in square brackets is the Fourier transform of the impact 
parameter space amplitude, U{x_i) — 1, where 

U{xA_) = Pexp (^-|^ A_{x+,x^)dx^^ (3.7) 

is the non-Abelian eikonal factor. The unit matrix 1 G SU{Nc), with Nc 
the number of colours, subtracts the field independent part, and the path 
ordering operator P sets the field at smallest x+ to the rightmost position. 

This formula was derived by many authors [28-30]. A derivation based 
on the summation of diagrams of the type shown in Fig. 9, is given in the 
Appendix of [6]. 

The amplitude of Eq. (3.6) is easily generalized to the case of a spinor 
quark or an energetic gluon. In the the high-energy limit, helicity-flip or 
polarization-flip contributions are suppressed. 



i2Tr5{ko-ko)T = 2Tr6iko-ko)2ko U{k'^ - k±) - i2Tr)^6^ik'^ - k±) 
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Figure 8: Scattering of a quark off the target colour field. 




Figure 9: Typical diagrammatic contribution to the cikonal amplitude, 
Eq. (3.6). Attachments of gluon lines with crosses correspond to vertices at 
which the classical external field appears. 

The eikonal approximation can be used for the calculation of the ampli- 
tude for qq pair production off a given target colour field [14] . Both diffrac- 
tive and inclusive cross sections are obtained from the same calculation, 
diffraction being defined by the requirement of colour neutrality of the pro- 
duced pair. 




Figure 10: Electroproduction of a qq pair off the target colour field. 

The process is illustrated in Fig. 10. The corresponding T matrix clement 
has three contributions, T = Tqq+Tq+Tq, where Tgg corresponds to both the 
quark and antiquark interacting with the field, while Tq and Tq correspond 
to only one of the partons interacting with the field. 

Let Vq{p',p) and Vq{k', k) be the effective vertices for an energetic quark 
and antiquark interacting with a soft gluonic field. For quarks with charge 
e, one has 

z27r(5(A;^, + fco - qQ)Tqq (3.8) 

= ze / ^^,,(y)^,(p^^');T^^(g)^^^#,fc')^r'(fc'), 
J (^TTj^ p — m —It — m 

where q = p + k hy momentum conservation, e(g) is the polarization vector 
of the incoming photon, and r', s' label the spins of the outgoing quarks. 
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The propagators in Eq. (3.8) can be treated in a high-energy approxi- 
mation. In a co-ordinate system where the photon momentum is directed 
along the z-axis, the large components are p+ and It is convenient to 
introduce, for each vector k, a vector k whose minus component satisfies 
the mass shell condition, k- = {k'j_ + m?)/k+, while the other components 
are identical to those of k. The propagators in Eq. (3.8) can be rewritten 
according to 



|/ — m — m? 2p_|_ 

and an analogous identity for the propagator with momentum k. In the high- 
energy limit, the term proportional to 7+ in Eqs. (3.9) can be dropped. 
After inserting Eqs. (3.9) into Eq. (3.8), the relation 

uAp')Vq{p' ,p)us{p) = 27rS{p'o-po)2po [u{p'±-p±) - {2TrfS''{p'^-p^)] S^s' , 

(3.10) 

which corresponds to Eq. (3.6), is applied. The vertex Vg{k,k') is treated 
analogously. Writing the loop integration as d^k = {\/2)dk+dk-(Pk± and 
using the approximation 5{lo) — 25{Ij^) for the energy (5-functions, the A;+ 
integration becomes trivial. The /c_ integral is done by closing the integra- 
tion contour in the upper or lower half of the complex k- plane. The result 
reads 

™ ie f ,9, a(l — a) , ^ ^ , ^ 

T,g = q+ j d^k^ ^2 + fc2_ ^ApyiqWik) (3.11) 



X 



U{p'± - P±) - i^T^fS^p'^ - p^)] \uHk± - fcl) - i27rf S\k'^ - k^) 



where = (1 — a)gf_|_, k'^ = ag+, A^^ = a(l — a)Q^ + m?. Thus, a and 
1— a characterize the fractions of the photon momentum carried by the two 
quarks, while {N'^+k\) measures the off-shellness of the partonic fluctuation 
before it hits the target. In the following we set m = 0. 

Adding Tq and Tq and introducing the fundamental function 

W^Ayi-) = U{xA_)U\x^ + y^)-l, (3.12) 

which encodes all the information about the external field, the complete 
amplitude can eventually be given in the form 



r = -^^+ / d'k^^2^''s'{p)^^q)vAk) l^e-'^'^^-W^Akl-kA, 

(3.13) 

where Wxj_ is the Fourier transform of Wxj_{y±) with respect to y±, and 
A± = k'j_+ p'j_ is the total transverse momentum of the final qq state. 
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Prom the above amplitude, the transverse and longitudinal virtual pho- 
ton cross sections are calculated in a straightforward manner using the ex- 
plicit formulae for Us'{p)^{q)vr'{k). Summing over all qq colour combina- 
tions, as appropriate for the inclusive DIS cross section, the following result 

is obtained, 



2e2Q2 
(27r)e 



dadk'1 
dadk'l ~ 2(27r)' 



{a{l-a)f j^^ j 



iV2 + fe2 



a) 



d'k 



k±W:,^{k'^-k^) 



N^ + k'] 



(3.14) 



(3.15) 



± 



The contraction of the colour indices of the two W matrices is implicit. 

An explicit calculation shows (see, e.g., [6]) that the leading contribu- 
tions at high Q2 are 

„2 



67r2Q2 



Jx± 



and 



ax 



(3.16) 



(3.17) 



^2j^J^iW,Ak'^-k^) 



N^+k"^: 



where < //^ < Q^. 

The resulting physical picture can be summarized as follows. For longi- 
tudinal photon polarization, the produced qq pair has small transverse size 
and shares the photon momentum approximately equally. Only the small 
distance structure of the target colour field, characterized by the quantity 
|9xVFa;j^(0)|2, is tested. For transverse photon polarization, an additional 
leading twist contribution comes from the region where a or 1 — a is small 
and k'^ ~ A^. In this region, the qq pair penetrating the target has large 
transverse size, and the large distance structure of the target colour field, 
characterized by the function Wx^iy±) at large y±, is tested. This phys- 
ical picture, known as the aligned jet model, was introduced in [12] on a 
qualitative level and was used more recently for a quantitative discussion of 
small- a; DIS in [13]. 

Diffractive cross sections are now derived by introducing a colour singlet 
projector into the underlying amplitude, i.e., by the substitution 



17 



in Eqs. (3.14) and (3.15). This change of the colour structure has crucial 
consequences. 

Firstly, the longitudinal cross section vanishes at leading twist (cf. 
Eq. (3.16)) since the derivative d±Wxj_{0) is in the Lie-algebra of SU{Nc), 
and therefore t'r d±Wxj_{0) = 0. 

Secondly, for the same reason the InQ^ term in the transverse cross sec- 
tion, given by Eq. (3.17), disappears. The whole cross section is dominated 
by the endpoints of the a integration, i.e., the aligned jet region. At leading 
order in the diffractive cross sections read 

af = (3.19) 



o2 



= —^-^ I da I dk'^ 



oo 



N^ + k"^, 



2 



.(3.20) 



The cutoff of the a integration, fi'^/Q'^, has been dropped since, due to the 
colour singlet projection, the integration is automatically dominated by the 
soft endpoint. 

In summary, the leading-twist cross section for small-a; DIS receives con- 
tributions from both small- and large-size qq pairs, the latter correspond- 
ing to aligned jet configurations. The requirement of colour neutrality in 
the final state suppresses the small-size contributions. Thus, leading twist 
diffraction is dominated by the production of pairs with large transverse size 
testing the non-perturbative large-distance structure of the target colour 
field. 



3.4- Diffractive parton distributions 

The basic theoretical ideas are due to Trentadue and Veneziano, who 
proposed to parametrize semi-inclusive hard processes in terms of 'fracture 
functions' [18], and to Berera and Soper, who defined similar quantities for 
the case of hard diffraction [19] and coined the term 'diffractive parton dis- 
tributions'. The following discussion is limited to the latter, more specialized 
framework. 

Recall first that a conventional parton distribution fi{y) describes the 
probability of finding, in a fast moving proton, a parton i with momentum 
fraction y. 

In short, diffractive parton distributions are conditional probabilities. 
A diffractive parton distribution dfP{y,^,t)/d$, dt describes the probability 
of finding, in a fast moving proton, a parton i with momentum fraction y, 
under the additional requirement that the proton remains intact while being 
scattered with invariant momentum transfer t and losing a small fraction 
^ of its longitudinal momentum. Thus, the corresponding 7*p cross section 
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can be written as [31] 



d^dt 



where a{x,Q^ ,y)'^*'^ is the total cross section for the scattering of a virtual 
photon characterized by x and and a parton of type i carrying a fraction 
y of the proton momentum. The above factorization formula holds in the 
limit DO with x, ^ and t fixed. 

At leading order and in the case of transverse photon polarization, only 
the quark distribution contributes. For one quark flavour with one unit of 
electric charge, the well-known partonic cross section reads 

giving rise to the diffractive cross section 

da{x, Q2, tr^P-P''' 2^e2 x df^{x, ^, t) 



d^dt Q2 ^^^i 



(3.23) 



where the factor 2 is introduced to account for the antiquark contribution. 

As in inclusive DIS, there are infrared divergences in the partonic cross 
sections and ultraviolet divergences in the parton distributions. They are 
conveniently renormalized with the MS prescription, which introduces the 
scale fi as a further argument. The distribution functions then read fi{x, /i^) 
and dfP{x,^,t,ii^)/d^dt in the inclusive and diffractive case respectively. 

Accordingly, Eq. (3.21) has to be read in the MS scheme, with a ji 
dependence appearing both in the parton distributions and in the partonic 
cross sections. The claim that Eq. (3.21) holds to all orders implies that 
these ^ dependences cancel, as is well known in the case of conventional 
parton distributions. Since the partonic cross sections are the same in both 
cases, the diffractive distributions obey the usual Altarelli-Parisi evolution 
equations, 

_J_ dfPix,C,t,^^) _ ^p..(,M&tl:lll^ (3 24) 

with the ordinary splitting functions Pij{x/y). 

Thus, for the analysis of diffractive DIS, it is essential to gain confidence 
in the validity of the factorization formula Eq. (3.21). Berera and Soper 
first pointed out [31] that such a factorization proof could be designed along 
the lines of related results for other QCD processes (see [32] for a review). 
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Proofs were given by Grazzini, Trentadue and Veneziano in the framework 
of a simple scalar model [33] and by Collins in full QCD [34]. 

The concept of diffractive parton distributions is more rigorous but less 
predictive than the older, widely used method of parton distributions of 
the pomeron [16, 17, 35]. This method originates in the observation that 
diffractive DIS can be understood as the soft high-energy scattering of the 
hadronic fluctuation of the photon and the target proton. Since a large 
sample of hadronic cross sections can be consistently parametrized using the 
concept of the Donnachie-Landshoff or soft pomeron [36] , it is only natural 
to assume that this concept can also be used in the present case. A more 
direct way of applying the concept of the soft pomeron to the phenomenon 
of hard diffraction was suggested by Ingelman and Schlein in the context of 
diffractive jet production in hadronic collisions [16]. Their idea of a partonic 
structure of the pomeron, which can be tested in hard processes, applies to 
the case of diffractive DIS as well [17]. Essentially, one assumes that the 
pomeron can, like a real hadron, be characterized by a parton distribution. 
This distribution factorizes from the pomeron trajectory and the pomeron- 
proton vertex, which are both obtained from the analysis of purely soft 
hadronic reactions. The above non-trivial assumptions, which have not been 
derived from QCD, are often referred to as 'Regge hypothesis' or 'Regge 
factorization'. 



3.5. Diffractive parton distributions in the semiclassical approach 

It will be shown how the semiclassical calculation can be factorized into 
hard and soft part and how a model for diffractive parton densities naturally 
arises from the soft part of the semiclassical calculation [37] . 

For simplicity, consider the fluctuation of the photon into a set of 
scalar partons which interact independently with the proton colour field 
(see Fig. 11). 




Figure 11: Hard diffractive process in the proton rest frame. The soft 
parton with momentum k is responsible for the leading twist behaviour of 
the cross section. 
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Assume furthermore that the transverse momenta p'^j-^^ij = l---n) are 
hard, i.e. 0{Q). A leading twist contribution to diffraction can arise only if 
the transverse momentum of the remaining parton is small, k'j^ ~ Aqcd- 

The standard cross section formula for the scattering off a static external 
field reads 



da 



-^|Tp27r%o-go)rf^^"+'^ where = A;' + EPm • (3-25) 

2qo 

All momenta arc given in the proton rest frame, T is the amplitude corre- 
sponding to Fig. 11, and dX^"'^^^ is the usual phase space element for n + 1 
particles. 

The colour rotation experienced by a each parton penetrating the exter- 
nal field is described by an eikonal factor. The resulting amplitude is given 

by 

i 2T,5{qQ -q'o)T = jTH 2TrS{k'o - ko) 2ko U{k'^ - k^ 

^^{jr^^^^iPm-Pim)MmU{p[^^^ , (3.26) 

where Tjj stands for the hard part of the diagram in Fig. 11. 

The integrations over the light-cone components P(j)+ can be performed 
using the energy (5-functions. The p(j)_-integrations are performed by pick- 
ing up the poles of the propagators Since the external field is smooth 
and Th is dominated by the hard scale, n — 1 of the n transverse momen- 
tum integrations can be performed trivially. This is not the case for the last 
integration which will necessarily be sensitive to the small off-shellness k'^. 
However, the n — 1 performed integrations ensure that the eikonal factors 
associated with the high-p^ partons are evaluated at the same transverse 
position. The resulting colour structure of the amplitude, after projection on 
a colour singlet final state, involves the trace of Wx^{y±) (cf. Eq. 3.12). It is 
intuitively clear that only two eikonal factors are present since the high-^j^ 
partons are close together in transverse space. They are colour rotated like 
a single parton. 

Under the further assumption that final state momenta of the high-p^ 
partons are not resolved on the soft scale, the following result is derived. 



da A;2|T^P " ' ' ^ ^ 



C JX_i_ 



tv[W^^{k'^-k^)] 
k± k"^ 



(3.27) 

In this expression the non-perturbative input encoded in the Fourier trans- 
form W is totally decoupled from the hard momenta that dominate Th- 
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Fig. 12. Diffractive DIS in the proton rest frame (left) and the Breit frame (right); 
asymmetric quark fluctuations correspond to diffractive quark scattering, asym- 
metric gluon fluctuations to diffractive boson-gluon fusion. 



The squared amplitude |Tffp in Eq. (3.27) can be expressed through 
the partonic cross section (T(a;, Q^, y). In Fig. 11 this corresponds to the 
interpretation of the line labelled by k as an incoming line for the hard 
process. The cross section differential in ^ takes the form 



Introducing the variable b = y/^ the diffractive parton density for scalars 
(integrated over t) reads 



£k^ tT[W^^{k'^-k^)] 



(27r)2 k'lb + kl{l-b) 



e\l-bj J (27r)4iV, 

(3.29) 

Analogous considerations with spinor and vector partons result in the same 
factorizing result, but with different distribution functions. 

The physical picture emerging from the above semiclassical calculation 
of diffractive parton distributions is illustrated in Fig. 12. In the Breit frame, 
leading order diffractive DIS is most naturally described by photon-quark 
scattering, with the quark coming from the diffractive parton distribution of 
the target hadron. This is illustrated on the r.h. side of Fig. 12a. Identifying 
the leading twist part of the qq pair production cross section (l.h. side of 
Fig. 12a) with the result of the conventional partonic calculation (r.h. side 
of Fig. 12a), the diffractive quark distribution of the target is expressed in 
terms of the color field dependent function given in Eq. (3.12). 
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Similarly, the cross section for the color singlet production of a qqg state 
(l.h. side of Fig. 12b) is identified with the boson-gluon fusion process based 
on the diffractive gluon distribution of the target (r.h. side of Fig. 12b). 
This allows for the calculation of the diffractive gluon distribution in terms 
of a function similar to Eq. (3.12) but with the U matrices in the adjoint 
representation. 

In the semiclassical approach, the cross sections for inclusive DIS arc 
obtained from the same calculations as in the diffractive case where, how- 
ever, the color singlet condition for the final state parton configuration is 
dropped. As a result, the qq production cross section (cf. the l.h. side of 
Fig. 12a) receives contributions from both the aligned jet and the high-p^ 
region. In the latter, the logarithmic dp\/p\ integration gives rise to a In 
term in the full cross section. 

In the leading order partonic analysis, the full cross section is described 
by photon-quark scattering. The gluon distribution is responsible for the 
scaling violations at small x, dF2{x^Q'^) / dhiQ"^ ~ xg{x, Q2). Thus, the 
semiclassical result for qq production, with its In contribution, is sufficient 
to calculate both the inclusive quark and the inclusive gluon distribution. 
The results are again expressed in terms of the function in Eq. (3.12) where 
now the color trace is taken after the two W matrices (corresponding to the 
amplitude and its complex conjugate) have been multiplied. 

To obtain explicit formulae for the above parton distributions, a model 
for the averaging over the color fields has to be introduced. Such models 
are provided, e.g., by the non-perturbativc stochastic vacuum [38] or by the 
perturbative small dipoles of [39]. The large hadron model [40], on which 
the following analysis is based, has the the advantage of justifying the Fock 
space expansion of the photon wave function while being intrinsically non- 
perturbative as far as the t channel colour exchange is concerned. 

In the case of a very large hadronic target [40] (see also [41]), even in the 
aligned jet region, the transverse separation of the qq pair remains small [42]. 
This is a result of the saturation of the dipolc cross section at smaller dipole 
size. Under the additional assumption that color fields in distant regions 
of the large target are uncorrelated, a simple Glauber-type exponentiation 
of the averaged local field strength results in explicit formulae for all the 
relevant parton distribution functions [15] (see [43] for a closely related 
analysis) . 

Thus, diffractive and inclusive quark and gluon distributions at some 
small scale Qq ^^e expressed in terms of only two parameters, the average 
color field strength and the total size of the large target hadron. The energy 
dependence arising from the large-momentum cutoff applied in the process 
of color field averaging can not be calculated from first principles. It is 
described by a In^ x ansatz, consistent with unitarity, which is universal for 
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Figure 13: The structure function F^^^\^, f3,Q'^) computed in the semi- 
classical approach with HI data from [46]. The open data points correspond 
to < 4 GeV^ and are not included in the fit. 



both the inclusive and diffractive structure function [44]. This introduces 
a further parameter, the unknown constant that comes with the logarithm 
(see [45] for a related but different way of introducing a phenomenological 
energy dependence in diffractive and inclusive DIS). 

A conventional leading order DGLAP analysis of data at small x and 
> Qo results in a good four parameter fit (Qo being the fourth param- 
eter) to both the inclusive and diffractive structure function (Figs. 13 and 
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Figure 14: The diffractive structure function /3, Q^) with data 

from ZEUS [47]. Open circles correspond to < 4 GeV^. The charm 
content is indicated as a dashed hne. 




14). Diffractive data with < 4 GeV^ is excluded from the fit since higher 
twist effects are expected to affect this region (cf. [24] for a phenomenolog- 
ical discussion of higher twist effects). As an illustration, the (3 dependence 
of F^^^^ at different values of is shown in Figs. 13 and 14 (see [15] for 
further plots, in particular of the inclusive structure function, and more 
details of the analysis) . 

Two important qualitative features of the approach should be empha- 
sized. First, the diffractive gluon distribution is much larger than the diffrac- 
tive quark distribution, a result reflected in the pattern of scaling violations 
of F^^^\ This feature is also present in the analysis of [39], where, in con- 
trast to the present approach, the target is modelled as a small color dipole. 
Second, the inclusive gluon distribution, calculated from qq pair production 
at high p± and determined by the small-distance structure of the color field, 
is large and leads to the dominance of inclusive over diffractive DIS. 



4. Vector meson production 

So far, inclusive diffraction, as parametrized, e.g., by the diffractive 
structure function ^2^, was at the centre of interest of this review. It was ar- 
gued that, for inclusive processes, the underlying colour singlet exchange is 
soft. In perturbative QCD, the simplest possibility of realizing colour singlet 
exchange is via two t channel gluons. In fact, the colour singlet exchange in 
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certain more exclusive diffractive processes is, with varying degree of rigour, 
argued to be governed by a hard scale. In such cases, two-gluon exchange 
dominates. 



4-1. Elastic meson production 

Elastic meson electroproduction is the first diffractive process that was 
claimed to be calculable in perturbative QCD [48,49]. It has since been 
considered by many authors, and a fair degree of understanding has been 
achieved as far as the perturbative calculability and the factorization of the 
relevant non-perturbative parton distributions and meson wave functions 
are concerned. 

To begin with, consider the electroproduction of a heavy qq bound state 
off a given classical colour field. The relevant amplitude is shown in Fig. 15. 
In the non-relativistic limit, the two outgoing quarks are on-shell, and each 
carries half of the J/^l^ momentum. Thus, the two quark propagators with 
momenta p' = k' = q' /2 and the J/V' vertex are replaced with the projection 
operator gj^j^^t' + m). Here 

Fgg is the electronic decay width of the J/V' particle, mj = 2m is its mass, 
and ej its polarization vector [50]. 




Figure 15: Leading order amplitude for the elastic production of a J/^* 

particle off an external colour field. 

Using the calculational technique of Sect. 3.3, the amplitude of Fig. 15 
can be expressed in terms of non-Abelian eikonal factors U and associ- 
ated with the two quarks. 

Considering both Q and m as hard scales, while U and U'^ are governed 
by the soft hadronic scale A, the integrand in the loop integration of Fig. 15 
can be expanded in powers of the soft momentum k±. The leading power 
of the amplitude is given by the first non-vanishing term. In the case of 
forward prodTiction, p'j_ = k'j_ = 0, the dependence on the external colour 
field takes the form [6] 



(fk±klti \uip'^-p±)U\k±-k'^) - {2'Kf5'^{p'^-pi_)5'^{k'^-ki_) 
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{i-KYd' / W.^(yx) _ . (4.2) 



3/X=0 



Now, the crucial observation is that precisely the same dependence on 
the external field is present in the amplitude for forward Compton scatter- 
ing shown in Fig. 16. In the case of longitudinal photon polarization, the 
transverse size of the qq pair is always small, and the target field enters only 
via the second derivative of W that appears in Eq. (4.2). 




Figure 16: The Compton scattering amplitude within the semiclassical 
approach. 

Thus, comparing with the parton model result for longitudinal photon 
scattering, this derivative can be identified in terms of the gluon distribution 
of the target proton [51], 

-dl^j trT^,^(2/^) =27rVx5(x). (4.3) 

Using this relation, the amplitudes for the forward production of trans- 
versely and longitudinally polarized J/^p mesons by transversely and lon- 
gitudinally polarized virtual photons are obtained. Under the additional 
assumption S> mj, the amplitudes for longitudinal and transverse po- 
larization read 

TL = -i647r^a,gjeixg{x))^ , Tt = '^Tl, (4.4) 

where ^/s is the ccntrc-of-mass energy of the 'y*p collision. This is Ryskin's 
celebrated result for elastic J/ip production [48]. 

It is not surprising that the gluon distribution of Eq. (4.3), calcu- 
lated according to Fig. 16, shows no scaling violations and only the triv- 
ial Bremsstrahlungs energy dependence ~ 1/x. The reason for this is the 
softness assumptions of the semiclassical calculation. Firstly, the eikonal 
approximation implies that all longitiidinal modes of the external field are 
much softer than the photon energy. Secondly, the reduction of the field 
dependence to a transverse derivative is only justified if the scales govern- 
ing the quark loop, i.e., in the case of Fig. 16 and and in the 
case of Fig. 15, are harder than the transverse structure of W. These two 
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approximations, evidently valid for a given soft field, are also justified for 
a dynamical target governed by QCD as long as only leading logarithmic 
accuracy in both 1/x and is required. Thus, a non-trivial dependence on 
1 /x and can be reintroduced into Eq. (4.4) via the measured gluon dis- 
tribution, keeping in mind that the result is only valid at double-leading-log 
accuracy. 

An essential extension of the above fundamental result is related to the 
treatment of the bound state produced. Brodsky et al. [49] showed that, 
at least for longitudinal photon polarization, a perturbative calculation is 
still possible in the case of light, non-perturbative bound states like the 
p meson. The calculation is based on the concept of the light-cone wave 
function of this meson (see, however, [53] for an alternative approach to 
light meson production). Referring the reader to [52] for a detailed review, 
a brief description of the main ideas involved is given below (cf. [54]). For 
this purpose, consider the generic diagram for the production of a light 
meson in a hard QCD process given in Fig. 17. 




Figure 17: Generic diagram for meson production in a hard process. 

Assume that, as shown in this figure, all diagrams can be cut across 
two quark-lines, the constituent quarks of the meson, in such a way as to 
separate the hard process H from the soft meson formation vertex V, which 
is defined to include the propagators. The amplitude can be written as 

T = j d^kTH{k)V{k)= j^dzTHiz)^ j dk-d^k±V{k) = dzTHiz)^iz), 

^ ° (4.5) 

where z = k-^/q'^, and the last equality is simply the definition of the light- 
cone wave function (f) of the meson. The two crucial observations leading to 
the first of these equalities are the approximate k- and k± independence of 
Tff and the restriction of the z integration to the interval from to 1. The 
first is the result of the hard scale that dominates Th, the second follows 
from the analytic structure of V. In QCD, the k± integration implicit in 
usually has an UV divergence due to gluon exchange between the quarks. 
Therefore, one should really read (f) = (j){z, ^^), where the cutoff is of the 
order of the hard scale that governs Th- At higher orders in a^, the hard 
amplitude Th develops a matching IR cutoff dependence. 

The discussion of vector meson production given so far was limited to 
the double-leading-log approximation as far as the colour singlet exchange 
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in the t channel is concerned. To go beyond this approximation, the concept 
of 'non-forward' or 'off-diagonal' parton distributions, introduced some time 
ago (see [55] and refs. therein) and discussed by Ji [56] and Radyushkin [57] 
in the present context, has to be used (see also [58]). 

Recall first that the semiclassical viewpoint of Figs. 15 and 16 is equiva- 
lent to two-gluon exchange as long as the transverse size of the energetic qq 
state is small. So far, the recoil of the target in longitudinal direction has 
been neglected. However, such a recoil is evidently required by the kinemat- 
ics. For what follows, it is convenient to use a frame where is the large 
component of the photon momentum. In Fig. 18, the exchanged gluons and 
the incoming and outgoing proton with momenta P and P' are labelled by 
their respective fractions of the plus component of P = (P + P')/2. If A is 
the momentum transferred by the proton, ^P_|_ = A+/2. The variable y is 
an integration variable in the gluon loop. 

Figure 18: Elastic J ftp production (three further diagrams with the gluons 
connected in different ways have to be added). The two gluon lines and the 
incoming and outgoing proton are labelled by their respective fractions of 
the plus component of P = (P -|- P')/2. 

The lower part of the diagram in Fig. 18 is a generalization of the con- 
ventional gluon distribution. It can be described by the non-forward gluon 
distribution 

HM^,i) = ^;^ J cix_e-2/P+-/2^p/|^t(o,:,_,o^)+/^P(0,0,Ox)/|P). 

The description of elastic meson production in terms of non-forward par- 
ton distributions is superior to the double-leading-log approach of [48, 49] 
since ag corrections to the hard amplitude, meson wave function and parton 
distribution function can, at least in principle, be systematically calculated. 
However, the direct relation to the measured conventional gluon distribution 
is lost. A new non-perturbative quantity, the non- forward gluon distribu- 
tion, is introduced, which has to be measured and the evolution of which has 
to be tested - a very complicated problem given the uncertainties of the ex- 
periment and of the meson wave functions involved (see [59] for possibilities 
of predicting the non- forward from the forward distribution functions). 
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Factorization 

Having discussed the leading order results for vector meson production, 
the next logical step is to ask whether the systematic calculation of higher 
order corrections is feasible. For this, it is necessary to understand the fac- 
torization properties of the hard amplitude and the two non-perturbative 
objects involved, i.e., the meson wave function and the non-forward gluon 
distribution. 

Factorization means that, to leading order in 1/Q, the amplitude can be 
written as 

T = £dzJ dyH{y,x/2)THiQ^y/x,z,iJ.'')(f>{z,fi^), (4.7) 

where Th is the hard scattering amplitude, (j) is the light-cone wave function 

of the vector meson produced, and H is the non-forward parton distribution 
of the proton. It could, for example, be the non-forward gluon distribution 
Hg of the last section. The variable x = xbj is the usual DIS Bjorken 
variable. 

In these notes, I will not describe the factorization proof for longitu- 
dinal vector meson production given in [60]. Instead, a particularly simple 
situation will be used to show how factorization works specifically at small 
X, from the point of view of the target rest frame [61]. For this purpose, 
consider a very energetic scalar photon that scatters off a hadronic target 
producing a scalar meson built from two scalar quarks (see Fig. 19). The 
quarks are coupled to the photon and the meson by point-like scalar ver- 
tices ie and iX, where e and A have dimension of mass. The coupling of the 
gluons to the scalar quarks is given by —ig {r/^ + r|^), where r and r' are the 
momenta of the directed quark lines, and g is the strong gauge coupling. 




a) b) c) 



Figure 19: The leading amplitude for a point-like meson vertex. 

Under quite general conditions [61], the gluon momenta satisfy the rela- 
tions i+, f+ <^q+,e^, £'_ < P_ and ~ £''^ £]_. Then the lower bubble 

in Fig. 19 effectively has the structure 



(4.8) 
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which is defined to include both gluon propagators and all colour factors. 

Assume that F restricts the gluon momentum to be soft, <C Q^. In 
the high-energy limit, it suffices to calculate 

where 

T^"" = Ti'^il, q) = T^" + T^"" + T^" (4.10) 

is the sum of the upper parts of the diagrams in Fig. 19. 

Note that, because of the symmetry of with respect to the two 
gluon lines, the amplitude T'^'^ of Eq. (4.10) is used instead of the properly- 
symmetrized upper amplitude 

irym(^,/,g) = \[T'^''ii,i',q)+T'''^{-i',-i,q)]. (4.11) 

The two exchanged gluons together form a colour singlet and so the sym- 
metrized amplitude satisfies the same Ward identity as for two photons, 

Ti^''^{£,e',q)e/, = 0. (4.12) 

Writing this equation in light-cone components and setting i± = £'_i, as 
appropriate for forward production, it follows that, for the relevant small 
values of ^_ , , i+ and i'^ , 

r,y^,++ ~ £l (4.13) 

in the limit — > 0. Here the fact that the tensor T^^, which is built from 
i', I and has no large minus components has been used. The i- integration 
makes this equation hold also for the original, unsymmetrized amplitude, 

y d^_r++ ~ l\ . (4.14) 

This is the crucial feature of the two-gluon amplitude that will simplify the 
calculation and lead to the factorizing result below. 

Consider first the contribution from diagram a) of Fig. 19 to the £- 
integral of T++ , which is required in Eq. (4.9) , 

Here = z{l — z)Q^, z = kj^/q-y. and the condition ^+ = 0, enforced by 
the (5-function in Eq. (4.8), has been anticipated. 
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Now / di-Th^^^ and / di-Tc^++ each carry no £± dependence. So, to 
ensure the vaUdity of Eq. (4.14), the sum of the three diagrams must be 

/d«-T,+ =4es\,/^.(l-.)Ar^3^, (4.16) 



where 



^f■■ 



7V2 + kl iV2 + {k± + ^i)2J (7V2 + A;2 )2 ■ 



(4.17) 



Note the l/Q'^ behaviour obtained after a cancellation of 1/Q^ contributions 
from the individual diagrams. This cancellation is closely related to the well- 
known effect of colour transparency. 

Introduce the k± dependent light-cone wave function of the meson 

Hz, kl) = / dk.dk^ (2;r)4fcV - ^^^^ " ^"^-^^^ 

The final result following from Eqs. (4.9) and (4.16) is a convolution of 
the production amplitude of two on-shell quarks and the light-cone wave 
function: 

^ = ''''' ( / ^^^^(^^)) I^^'I ''"'^WTWf'^^'^ • ^'-''^ 

This corresponds to the O(^j^) term in the Taylor expansion of the contri- 
bution from Fig. 19a, given in Eq. (4.15). 




Figure 20: Diagram for meson production with the vertex modelled by 
scalar particle exchange. 

At leading order, factorization of the meson wave function was trivial 
since the point-like quark-quark-meson vertex V{k'^, {q' — k)"^) = iX was 
necessarily located to the right of the all other interactions. To see how 
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factorization comes about in the simplest non-trivial situation, consider the 
vector meson vertex 

which corresponds to the triangle on the r.h. side of Fig. 20. Here, the dashed 
line denotes a colourless scalar coupled to the scalar quarks with coupling 
strength A'. 

The diagram of Fig. 20 by itself gives no consistent description of meson 
production since it lacks gauge invariance. This problem is not cured by just 
adding the two diagrams 19b) and c) with the blob replaced by the vertex 
V. It is necessary to include all the diagrams shown in Fig. 21. 




Figure 21: The remaining diagrams contributing to meson production 
within the above simple model for the meson wave function. 

The same gauge invariance arguments that lead to Eq. (4.14) apply to 
the sum of all the diagrams in Figs. 20 and 21. Therefore, the complete result 
for which is now defined by the sum of the upper parts of all these 

diagrams, can be obtained by extracting the term at leading order in the 
energy and Q^. Such a term, with a power behaviour ~ i'j^/Q^, is obtained 
from the diagram in Fig. 20 (replace i\ in Eq. (4.15) with the vertex V 
of Eq. (4.20)) by expanding around i± = 0. It can be demonstrated that 
none of the other diagrams gives rise to such a leading-order £\ contribution 
(see [61] for more details). 
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The complete answer is given by the term from the Taylor expansion 
of Eq. (4.15). The amplitude M is precisely the one of Eqs. (4.19) and 
(4.18), with i\ substituted by V of Eq. (4.20). The correctness of this simple 
factorizing result has also been checked by explicitly calculating all diagrams 
of Fig. 21. 

The above simple model calculation can be summarized as follows. The 
complete result contains leading contributions from diagrams that cannot 
be factorized into quark-pair production and meson formation. However, 
the answer to the calculation can be anticipated by looking only at one 
particular factorizing diagram. The reason for this simplification is gauge 
invariance. In the dominant region, where the transverse momentum £_l of 
the two i-channel gluons is small, gauge invariance requires the complete 
quark part of the amplitude to be proportional to l^. The leading i\_ de- 
pendence comes exclusively from one diagram. Thus, the complete answer 
can be obtained from this particular diagram, which has the property of 
factorizing explicitly if the two quark lines are cut. The resulting amplitude 
can be written in a factorized form. 

Note finally that, although QCD predicts factorization for longitudinal 
vector meson production and a l/Q^ suppression of the transverse cross 
section, this behaviour is not visible in the data [62] . A recent calculation [63] 
explaining the data is, in my opinion, in conflict with QCD expectations 
and, in particular, with the implications of [61], since it employs a non-local 
vector meson vertex in a loop calculation. It is important to clarify this 
situation. 

5. Summary 

The phenomenon of diffraction in DIS can be understood in very simple 
terms if the process is viewed in the rest frame of the target proton. The 
energetic virtual photon develops a partonic fluctuation which then scatters 
off the target. A certain fraction of the total DIS cross section is due to 
photon fluctuations with large transverse size. This fraction is not power 
suppressed in the high-Q^ limit. Since, as one expects from the experience 
with hadron-hadron collisions, part of these large size fluctuations scatter 
quasi-elastically off the proton, a leading twist diffractive cross section is 
obtained. 

One possibility to describe the interaction of the photon fluctuation with 
the target is two-gluon exchange. In certain more exclusive processes, such 
as longitudinally polarized vector meson production, the transverse size of 
the relevant photon fluctuations is always small and two-gluon calculations 
can be rigorously justified. However, for most of the diffractive cross section 
no suppression of multiple gluon exchange can be derived. 
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This problem is addressed in the semiclassical approach, where the tar- 
get is described by a superposition of soft colour fields and gluon exchange 
is resummed to all orders in an eikonal approximation. Diffraction arises 
whenever the partonic fluctuation of the photon remains in a colour sin- 
glet. The application of the semiclassical approach to experimental results 
is particularly simple if the approach is used to derive both diffractive and 
inclusive parton distributions at some small input scale. In this case, the 
analysis of all higher-Q^ data proceeds with standard perturbative meth- 
ods. Different models for the underlying colour fields can be compared to 
diffractive and inclusive structure function data in a very direct way. 
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